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How does the classical phase space structure for a composite system relate to the entanglement 
characteristics of the corresponding quantum system? We demonstrate how the entanglement in 
nonlinear bipartite systems can be associated with a fixed point bifurcation in the classical dynamics. 
Using the example of coupled giant spins we show that when a fixed point undergoes a supercritical 
pitchfork bifurcation, the corresponding quantum state - the ground state - achieves its maximum 
amount of entanglement near the critical point. We conjecture that this will be a generic feature of 
systems whose classical limit exhibits such a bifurcation. 



I. INTRODUCTION 

With the advent of quantum information theory, en- 
tanglement is now regarded as a physical resource that 
can be utilized to perform numerous quantum compu- 
tational and communication tasks This has in turn 
led to the study of the entanglement characteristics of 
various systems, and in turn, how these characteristics 
relate to more fully understood properties of the system. 
Such studies are two-fold beneficial - further elucidating 
the nature of entanglement as well as providing a new 
approach to the study of complex, quantum many-body 
systems. 

One area where such an approach has had some suc- 
cess is in the study of quantum phase transitions (QPTs) - 
qualitative changes in the ground state of a multi-partite 
system induced by the variation of some external param- 
eter 0. There have been many recent studies relating 
entanglement and QPTs (see U II HIE II 110) ■ 

Generally it has been found that in infinite systems 
that undergo a quantum phase transition at a critical 
parameter value, A — Ac, the entanglement as a function 
of A is a maximum at Ac- Several examples include: (i) 
The single site entanglement and the next nearest neigh- 
bor concurrence of the transverse Ising chain |1,^0 (al- 
though the nearest neighbor concurrence does not have 
its maximum value at A = Ac, its first derivative with 
respect to A does Q), (ii) the entropy of entanglement of 
half of a XXZ spin chain in a magnetic field ^ , and (iii) 
the entropy of entanglement of a single qubit with a bath 
of oscillators (the spin-boson model) 0. Such systems 
demonstrate a correspondence between quantum critical 
phenomena and entanglement. 

A QPT corresponds to a qualitative change in the 
ground state as a system parameter is varied. In the clas- 
sical regime, minimum energy coordinates correspond to 
elliptic (stable) fixed points. As a parameter in the sys- 
tem is varied, fixed points may undergo bifurcation jlOj - 
a loss of stability, the emergence of new fixed points - at 
some critical value of the parameter. This corresponds 
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to a qualitative change in the phase space structure of 
the system. 

In this article we consider the ground state of a system 
whose classical limit exhibits a bifurcation where a single 
elliptic fixed point loses its stability while two new, el- 
liptic points emerge - a so-called supercritical pitchfork 
bifurcation. Elliptic fixed points can be associated with 
the ground state of the quantized system. Subsequently 
we expect to see some signature of the classical bifur- 
cation in the quantum ground state, around the critical 
point. We argue that this signature is a peak in the en- 
tanglement with respect to the bifurcation parameter. 

Schneider and Milburn alluded to such a correspon- 
dence in their work on the Dicke model 11]. It was shown 
that the entanglement in the steady state of this system 
is a maximum for the parameter value corresponding to a 
bifurcation of the fixed points in the corresponding clas- 
sical dynamics. It was conjectured that the loss of sta- 
bility of a classical fixed point due to such a bifurcation 
will generically be associated with entanglement in the 
steady state of the full quantum system. We show that 
it is specifically the pitchfork nature of this bifurcation 
that is responsible for the peak in the ground state en- 
tanglement. 

To demonstrate this, we use the example of coupled gi- 
ant spins. This system is motivated by a proposed phys- 
ical implementation for quantum computation |l2|. In 
this proposal, qubits are realized by magnetic clusters 
- nanometer scale molecular clusters that have all the 
attributes of mesoscopic systems such as angular mo- 
mentum and magnetic moment. Two qubit gates are 
constructed by the coupling of the clusters via supercon- 
ducting quantum interference devices (or SQUIDs), as 
shown in figure ^ 

The strength of the coupling is dependent upon the 
super-current induced in the loop by one spin and the 
field this produces at the other site 12]. This field in- 
duces a rotation about the 'z'-axis of one cluster with 
frequency proportional to the z-component of the angu- 
lar momentum of the other cluster. In this way, we can 
imagine the system as a set of 'spinning tops', coupled 
via a nonlinear interaction, described by Hint = Jz® Jz- 

The article is organized as follows: In Sec. ^l^e intro- 
duce the coupled tops model. We begin with a review of 
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FIG. 1: (Color online) A schematic diagram of the cou- 
pled qubit realization of Ref. 12]. The magnetic clusters 
(the qubits) are coupled to superconducting loops of micro- 
SQUIDS and arranged in a ID lattice. Josephson junction 
switches are used in the coupling circuits, as shown. 



the classical dynamics based on Skellett and Holmes [I^ , 
focussing on the bifurcation analysis. This is followed by 
a calculation of the ground state entanglement for the 
quantum system, as the bifurcation parameter is varied. 
Here the link between the ground state entanglement and 
the bifurcation is first established. Using the Husimi dis- 
tribution (or Q-function) to represent the ground state 
in phase space, we demonstrate how the structure of the 
ground state is related to the classical fixed points. The 
aim of Sec. IHIl is to explain the observed behaviour and 
argue its generality. We finish with a discussion of our 
results and future directions. 



II. THE COUPLED TOPS MODEL 

The simplest Hamiltonian describing the coupled tops 
is a generalization of the N = 2 case of the transverse 
field quantum Ising model, i.e. 

H = (jjjx® i + ® Jx + -Jz® Jz (1) 

i 

where the angular momentum operators Ja satisfy the 
SU{2) commutation relations [JxT-Jy] — iJz (and cyclic 
permutations) . Setting /i = ^ , results in a one-parameter 
Hamiltonian 



H — Jxl + Jx2 + — JzlJz2 

J 



(2) 



where we make use of the notation Jai = Ja ® I and 
Ja2 = I ® Ja, a, = x,y,z, such that the subscript 1 (2) 
refers to subsystem 1 (2). 

The square of the total angular momentum of the in- 
dividual tops, 



satisfy 



Ji,H — J2,H — 0. 



(3) 
(4) 



so are constants of the motion. Here the tops are identi- 
cal, such that Jf = J| ~ + This allows the system 
to be represented in the basis of tensor products of the Jz 
eigenstates, |j, to)® \ j, n) = |to, n), where —j < m,n < j. 
Note that the coupling term in the Hamiltonian is scaled 
with j, to allow the classical limit to be taken. For j = ^, 
the Hamiltonian Q is analogous to the quantum Ising 
model for 2 spins, studied in Ref. 

Interestingly, the square of the total angular momen- 
tum of the system. 



J2^ j2 + j2 + 2Ji-J2, 



(5) 



is not a constant of the motion. The reason is that an 
external control is required to couple the tops. In our 
motivating example this is the circuit that inductively 
couples the tops. This is similar to the situation for two 
interacting qubits. A general two qubit gate, such as a 
controlled-NOT gate does not conserve total angular mo- 
mentum either (though the gate can fix the singlet and 
triplet subspaces). Again this is due to external interac- 
tions that control the gate. 

The classical analogue of this system has been rigor- 
ously studied by Skellett and Holmes fTS]. We now de- 
rive the semiclassical limit of Hamiltonian (|2J), showing 
its correspondence to the model of Ref. J.3J, and review 
the key points of this analysis relevant to our work. 



A. Classical Description 

The semiclassical limit of the coupled tops system cor- 
responds to the limit of j ^ oo. To obtain the semi- 
classical model, we express the classical coordinates as 
Laa = {Jaa) I j- In the limit, this allows the factorization 
of all moments, i.e. {JxiJzi) / P ~^ LxiL^i (for details 
see Appendix A). The classical equations of motion are 
obtained from the Heisenberg operator equations of mo- 
tion by taking expectation values and applying the fac- 
torization rule above. The corresponding semiclassical 
Hamiltonian is 



E — Lxi + Lx2 + tJ'LziLz2 



(6) 
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where E = {H)/j with the spherical constraint, L^^ -I- 

zi I I 

The analysis of the corresponding classical system [13| 
has shown that the non-linearity of the interaction term 
leads to chaotic motion for given parameter ranges and 
initial coirditioirs. More relevant in this context is the ex- 
istence of a supercritical pitchfork bifurcation, at a criti- 
cal value of the coupling parameter. From Ref.^3| for the 
semiclassical system, the critical value is /ic = 1. Below 
this critical value the dynamics of the system is predomi- 
nantly regular while above the phase space is mixed, with 
extensive regions of chaotic motion. 

The fixed points of the system are found by setting the 
equations of motion to zero, to determine the coordinates 
where the phase space flow is zero. From Skellett and 
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solely in the directions, denoted (IT), (Ti)j (it) and (|,|). 



FIG. 2: (Color online) The pitchfork bifurcation at /ic = 1 for 
Z/zi. The elliptic (stable) fixed points are centres, indicating 
that close to the fixed point, the motion is periodic, while the 
hyperbolic (unstable) point is a saddle. This diagram is the 
same for both values of L^i = = ±1. 



Holmes there are four solutions which exist for all 
values of the coupling parameter /x, given by 



L. 



■y2 



0. (7) 



At the critical value, = 1 the two fixed points at L^x = 
= 1 and = L^i = — 1 bifurcate, resulting in the 
emergence of a further jour fixed points, located at 



I^x\ — Lx2 — — , Lzi — — ±\/ 1 ^1 

[I V M 

— Ly2 — 0, 



Lxl — Lx2 — , Lzl — —Lz2 — n^ 

fl V M 

Lyl — Ly2 =, 



(8) 
(9) 



which exist for all /i > 1 . The stability of the fixed points 
is determined by analysis of the eigenvalues of the lin- 
earized matrix about each fixed point [lol |. 

In Ref. , it was shown that the two fixed points (O 
with Lxl = —Lx2 are unstable for all values of /i. The 
points, Lxl = Lx2 = ±1 are stable for /i < 1, becom- 
ing unstable at ^.c- The emergent fixed points are all 
stable. This implies that the bifurcations occurring at 
Lxl = Lx2 = il are supercritical pitchfork bifurcations, 
as illustrated in figure [21 

Since the total angular momenta of the two tops are 
conserved, their dynamics of each are constrained to the 
unit sphere, in angular momentum space. It is possible 
to reformulate the dynamics in terms of spherical polar 
coordinates - the polar angle from the positive Lzi axis, 
< < TT, and the azimuthal angle in the Lxi — Lyi 
plane (from the positive Lxi axis), < 0^ < 2ti. These 
coordinates give the angular momentum components via 



Lxi = sin 9i cos < 



sm Vi sm ( 
cos 9i . 



(10) 



For all fixed points Lyi — Ly2 — 0, which in spherical 
polar coordinates, corresponds to (f>i,(t>2 = or tt. Thus, 
we can view the fixed points as lying on the unit circle in 
the Lxi ~ Lzi planes, characterized by the polar angles, 
9i. For /i below the critical coupling there are two, stable 
fixed points, both of which lie at the 'equator' of these 
unit circles {9i = ^) at Lxi = Lx2 = 1 (0i = 4'2 = 0) and 
at Lxl = Lx2 = —1 {(^i = 4>2 = tt). Following the nota- 
tion used by Skellett and Holmes [13 , we denote the two 
states by (—>—>) and (^^) respectively, corresponding 
to the direction of the angular momentum vector. 

For ^ greater than the critical value, there are four 
stable fixed points, whose positions are shown in figure 
O The points labeled A and B correspond to Eq. © 
and points C and D to Eq. ©. Clearly as ^ ^ 00, 

—> (in figure O and, in the pictorial (arrow) notation 
of above, we have the four fixed points at the 'poles' of 
the spheres, in the four combinations (tt), (Ti)j (it) and 
(II), all of which are stable. 

The semiclassical analysis we have presented here is 
just a brief summary of those aspects most relevant for 
this paper. For an in-depth analysis of the classical dy- 
namics of this system we again refer the reader to Skellett 
and Holmes p!3 |. 

The bifurcating fixed point at Lxi — Lx2 = — 1 (for 
/i < /ic) is the minimum energy point and so corresponds 
to the quantum ground state. We now consider the quan- 
tum regime, and the entanglement between the spins as 
a function of the coupling strength, /u. 



B. Ground State Entanglement 

The ground state is computed by direct numerical di- 
agonalization of the Hamiltonian |(5J). The entanglement 
measure we employ is the entanglement of formation, 
which, for pure bipartite states, is equivalent to the en- 
tropy of entanglement 1, 15] 



S{Pi) = -Tr [pilogpi) 



ill) 



where pi = Tr^ (p) is the reduced density operator and 
the log is to base 2. 
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FIG. 4: (Color online) Variation in the entropy of entangle- 
ment of the ground state for with respect to the coupling 
strength n and the total subsystem angular momentum, j. 
Note that as the system becomes more classical as j increases 
that the peak in the entanglement versus /i becomes more 
evident. 
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FIG. 5: (Color online) The variation in the quantum critical 
parameter, ^iqc, the value of pi for maximum ground state en- 
tanglement, with total angular momentum j. In the classical 
limit, J — > CX3, approaches the classical critical parameter, 
where the bifurcation occurs. 



The case of j = 1/2, i.e. a two-site transverse field Ising 
model, was considered by Gunlycke et al. It was 

shown in this case that the ground state entanglement 
is zero for zero coupling (/i = 0) then increases as the 
coupling increases, to asymptote to the maximal value of 
1 as /i — > oo. 

For j > 1/2, the entanglement of the ground state 
with respect to the interaction strength, ji takes on a 
new characteristic (see figure 

For /i = 0, the ground state is simply a tensor prod- 
uct of the minimal J^i weight states, | — j,—j)x^ which 
is separable. As /i — + cxd, the ground state approaches 
the superposition, (|j, —j) + \ — /\/2, so the entan- 
glement still asymptotes to 1. However for j > 1/2 the 
entanglement now peaks at a finite value of /i. The height 
of this peak grows with the value of j and will approach 
infinity in the limiting case. 

We let fiqc denote the quantum critical parameter, de- 
fined as the coupling value at which the ground state 
entanglement is maximum. Figure El shows the limiting 
behaviour of fXqc- For larger j, the maximum entangle- 
ment occurs near the bifurcation point and in the limit 
j —^ oo, the quantum critical point approaches the clas- 
sical bifurcation point, 

fJ-qc ^J'C■ (12) 

In the semiclassical limit, the quantum critical point cor- 
responding to maximum ground state entanglement is 
the classical bifurcation point. To understand how the 
fixed point structure manifests in the quantum regime, 
we need to consider the structure of the quantum state 
in phase space. 



C. Ground State in Phase Space 

Fundamental to any comparison of classical and quan- 
tum dynamics is some notion of the quantum analogue of 
a classical joint phase-space probability distribution. We 
choose the Husimi or Q-function as the appropriate 
quantum analogue of the classical phase-space density, 
following ref. 's [ll El IH IS HI El . 

For systems described in spherical phase-space, Ap- 
pleby 23] demonstrated that the positive operator valued 
measurement (POVM) for optimal simultaneous mea- 
surements of angular momentum components is given by 
£^(2:) = |2:)(2|, where \z) are the SU{2), coherent states 

\z)^il + \z\T^e^J^\j,-j) (13) 

with |j, ~j) the lowest weight eigenstate of Jz and z, the 
stereographic projection of the sphere onto the plane, 

z(6»,<?!)) = e-'-^tan^. (14) 

The probability distribution for measurements defined by 
this POVM is the angular momentum representation Q- 
function, which for pure state p — is 

Q^(z)=rr(p£;(z)) =|(z|^)|^ (15) 

For our bipartite system we define the two-body Q- 
function in terms of the coherent states 

\z) = \zi)^\z2). (16) 
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FIG. 6: (Color online) The (j>i = (f>2 = tt cross-sections of the ground state Q-function for a coupling strength /i equal to 
(a) 0, (b) 0.7, (c) 1, (d) 1.184, (e) 1.25 and (f) 1.55, for j = 14. The x and y axes correspond to 9i and 02 respectively. 
The entanglement of formation is for (a), approximately 1 for (f), and maximum for (d), when the distribution is the most 
delocalized. We see that as the fixed point bifurcates, so does the ground state Q-distribution, from the single-peak to the 
twin-peaked structure. 



Since the Q-function for the coupled tops is in 4 phase- 
space dimensions, we may only display cross-sections 
graphically. We have calculated the (j)! = 'p2 — cross- 
section of the ground state Q-function for various values 
of the coupling strength, /i. This cross-section was cho- 
sen since this is the plane in which the bifurcating fixed 
points are located. 

As evidenced in figure El the structure of the ground 
state in phase space is intimately related to the fixed 
point structure and bifurcation. For small fi, the Q- 
distribution of the ground state is localized around the 
fixed point at di = 62 = n/2 which is elliptic for 
fi < j-ic- For /i much larger, beyond the bifurcation, the 
Q-distribution is twin-peaked, localized around the two 
emergent elliptic fixed points. Put simply, in these two 
extremes, the distribution is localized around the ellip- 
tic fixed points. In between, the distribution is spread 
between the three fixed points and it is this region in 
which the degree of entanglement is greatest. In fact, we 
will argue that the entanglement is maximum when the 
distribution is at its most delocalized psj . 

Under this assumption, the maximum ground state en- 



tanglement will occur when /i is above, but close to the 
critical point. This is when the Q-distribution is at its 
most delocalized, smeared between the three fixed points. 
As j becomes very large, and the Heisenberg uncertainty 
limit allows for a higher degree of localization of the dis- 
tributions, this quantum critical point approaches the 
classical bifurcation point. In this semiclassical regime, 
the greatest delocalization occurs closer and closer to the 
bifurcation, as — > 00. 



For this simple coupled tops model we have demon- 
strated that the entanglement characteristics of the 
ground state can be associated with the supercritical 
pitchfork bifurcation. To validate and generalize this 
result, we need to consider why the ground state Q- 
distribution is related to the fixed point structure, and 
how this structure corresponds to the degree of entangle- 
ment. 
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III. THE GROUND STATE Q-FUNCTION AND 
ENTANGLEMENT 



There has been considerable effort in relating quan- 
tum eigenstates to classical phase space structures (for 
example see |23, H^). The Einstein- Br illouin- Keller 
(EBK) approximation '2911 provides for a semiclassical 
quantization of classically integrable systems, whereby 
eigenstates are identified with closed loops around in- 
variant tori in phase space. For non-integrable systems, 
the Gutzwiller trace formula allows for semiclassical ap- 
proximations of energy spectra, but not the structure of 
the corresponding eigenstates |2^. 

For systems exhibiting mixed phase space, while no 
general quantization procedure exists, there is an under- 
standing that eigenstates can be separated into regular 
and irregular/chaotic |3Cll | (and in-between or hierarchical 
l^s'l) groups. Regular states are supported by classical 
tori obeying EBK quantization, whereas chaotic states 
are associated with chaotic phase space regions. In the 
Husimi representation, regular states are seen as local- 
ized on regular trajectories, whereas chaotic states are 
somewhat evenly distributed over the chaotic region of 
phase space [3ll |. 

Classically, the minimum energy trajectories corre- 
spond to elliptic (stable) fixed points. Thus, it would 
seem reasonable to infer that the quantum ground state 
would be localized around such points. 

The aim of this section is to demonstrate that for 
Hamiltonians whose classical analogue exhibits (mini- 
mal energy) elliptic fixed points, the Q-function of the 
ground state is peaked on the phase space coordinates 
corresponding to such fixed points. This statement has 
non-trivial consequences when the classical Hamiltonian 
displays a bifurcation of fixed points, which, in the quan- 
tum regime, will correspond to a qualitative change in the 
structure of the ground state phase-space distribution. 
Subsequently, we relate the structure of the Q-function 
to entanglement properties. 



A. Ground state Q-function and fixed points 

In reformulating quantum mechanics in phase space, 
the coherent states provide a natural phase space struc- 
ture for a given quantum system, as well as useful distri- 
butions based in the coherent state representation |32l |. 
two of which, we briefly discuss now. 

An operator O can be expressed in the diagonal form 



O = I Op{z)\z){z\dn{z) 



(17) 



called the P representation, where \z) (in complex vari- 
able z) are the coherent states and dii{z) is the corre- 
sponding measure. The density operator is given by 



p= / P{z)\z){z\dpi{z). 



(18) 



The Q-representation of an operator is defined as 

6 Oq{z) ^ {z\d\z) (19) 

and the density operator in the representation is denoted 
by Q(z). The statistical average of an operator in the 
two representations is given by 



(O) = Tr(pO) = / Qiz)Op{z)d^i{z) 
Piz)OQ{z)d^i{z). 



To determine the ground state of a given Hamiltonian, 
H, in this phase space formulation we seek the distribu- 
tion that minimizes the average energy, i.e. determine 
the distribution Q{z) that minimizes 



{H) = / Q{z)Hp{z)dpi{z) 



(20) 



for a given Hp{z). 

Upper and lower bounds on the ground state energy, 
Eq, are determined by minimizing the Q and P represen- 
tations of the Hamiltonian |32|. 



min_ffQ(z) > Eq > min Hp{z). 



(21) 



In the semi-classical limit (i.e ft — > for Heisenberg- 
Weyl, J — > oo for SU{2) coherent states) of the appropri- 
ately scaled Hamiltonian, both the Q and P representa- 
tions converge to the classical Hamiltonian, H{z). So the 
problem of finding the ground state in the phase space 
representation as we approach the semiclassical limit re- 
duces to finding the distribution Q{z), minimizing the 
functional, 



T[Q{z)] ^ / d^l{z) Q{z)n{z) 



(22) 



where 'H{z) is the semiclassical Hamiltonian. This distri- 
bution must satisfy 



dii{z) Q{z)^l, 0<Q(z)<l 



(23) 



as well as the uncertainty principle, which excludes the 
possibility of delta functions. 

Let us consider a suitably scaled Hamiltonian with 
semiclassical limit, Ti.{z), which exhibits a global mini- 
mum at z — zq. This requires that 



m 

dz 



and 



dzdz 



m 

dz 



dz"^ 



= 



{dzf 



> 0. 



(24) 



(25) 
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This corresponds to an elliptic fixed point in the phase 
space portrait of the classical dynamics. Taylor expand- 
ing H{z) about the fixed point gives 



H{z,z) = 7i(zo, Zq) + |z - zoP 



dzdz 



1/ X2 



h.o.t. 



which can be written as 



where 



H(zo, Zq 


) + [^ 


- Zq Z- 
















dzdz 








Zo 


^0 












dzdz 


Zo ^ ' 


Zo - 



-2n 



(26) 



such that det(D) < 0. In other words, to second order 
around the elliptic fixed point, 7i(2:, z) is parabolic, with 
minimum at zq. 

Now consider the minimization problem of Eq. H22|l . If 
the system under investigation has only a solitary global 
minima, since Q{z) takes only positive values, the distri- 
bution that minimizes Eq. H22|) will be as sharply peaked 
as possible on the coordinates of the minima. The ground 
state Q-function will be peaked on the coordinates of the 
fixed point. 

What if there are more than one (degenerate) global 
minima? In this case, consider expansions around each 
of the fixed points Zfc following (j^ . We are interested in 
the ground state in a semiclassical regime where the Q- 
functions can be highly localized. We thus approximate 
7i as a piecewise sum of the expansions around each of 
the degenerate fixed points. The functional H22|) is thus 
written as the summation 



1. JTlk 



d^i{z) Qk{z)nk{z) (27) 



where Ti.k{z) is the expansion around the fixed point Zk 
and Qk{z) is the distribution defined in the phase-space 
region TZk- Minimizing each term in this summation gives 
in turn a piecewise defined distribution, sharply peaked 
around the coordinates of each of the fixed points. 

The above statements are by no means rigorous. We 
are considering the problem in a semiclassical regime, 
which allows us to take simple approximations for the 
Hamiltonian, and simplifies the structure of the corre- 
sponding ground state Q-functions. However, as we take 
the semiclassical limit (ft — ^ 0, j — > oo), the Q- function 
will merge smoothly from a true quantum ground state 
distribution to the more semiclassical structure described 
above. 



As we move in the other direction, from the semiclas- 
sical to the quantum, we expect to observe remnants of 
the semiclassical structure. The Q-function, while per- 
haps not sharply peaked, will be concentrated around the 
fixed point coordinates. 

These observations lead us to make the following con- 
jecture: 

Conjecture 1 (Fixed point correspondence) Let 

H he a Hamiltonian whose classical analogue is defined 
as 'H{x,X*) = (xl^lx); where {\x)}, is the set of 
coherent states corresponding to the topology of the 
phase-space of Ti.{x,X*)- Then the Husimi distribution 
of the ground state of H will be concentrated around 
the phase-space coordinates corresponding to the fixed 
point (s) ofTi. 

While we have not provided a rigorous proof of this state- 
ment, we have argued above that such a notion is indeed 
plausible. 

This conjecture by no means provides detailed informa- 
tion about the ground state, but instead provides a guide- 
line for the phase-space structure of the ground state. As 
in Ref. [s^l , where the use of the Q-representation of the 
Hamiltonian as a guide to studying ground state phase 
transitions was advocated, this guideline become most 
useful when the fixed points vary according to some crit- 
ical parameter i.e. bifurcation. This can provide for a 
variational approach to approximating the ground state 
via coherent states (see Ref. H^l)- 

For our purposes, this link between fixed points and 
the phase-space structure can provide a qualitative un- 
derstanding of the effect of the bifurcation on the ground 
state entanglement. The connection to entanglement is 
made via the work of Sugita 25], which we now discuss. 



B. Entanglement and the Q-function 

When dealing with Q-functions for composite systems, 
there are several ways to generalize the idea of coherent 
states to many- body systems. In the next section, we 
follow Sugita (23 by constructing coherent states based 
on the single-particle transformation group. This means 
the coherent states are independent-particle states, and 
hence separable, allowing the entanglement to be related 
to the structure of the Q-function. 

Sugita [25 constructs the many-body coherent states 
based on the single-particle transformation group, fol- 
lowing the group-theoretical construction of Perelomov 
[2^ 1. As an illustrative example, consider a system of 
n qubits (basis states | t),| i))- The local unitary (or 
single-particle) transformation group is S'[/(2)®". Co- 
herent states are generated by applying this group to the 
lowest weight state | |)®". In this way, we generate all 
(pure) separable states. From this definition, coherent 
states are thus equivalent to separable states. We note 
that our arguments in the above section easily apply to 
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the multi-party case when the coherent states are gener- 
ahzcd in this way. 

From this definition, Sugita argues that a separable 
state is represented by a localized wave packet in phase 
space. Since coherent states are the most localized states 
in the Husimi representation, it is argued that delocal- 
ization of the Husimi distribution implies correlation - 
hence entanglement - between the particles. This delo- 
calization can be measured by the Renyi-Wehrl entropy, 
which represents the effective volume occupied by the 
Husimi distribution. 

The majority of Ref. [2^ is devoted to constructing 
explicit formulas for the Renyi-Wehrl entropy in terms of 
the moments. The calculation of the moments is simpli- 
fied via a group-theoretical construction. When applied 
to a system of two qubits, it's shown that the concur- 
rence iTiay be expressed in terms of the second mo- 
ment of the Q-function. More strikingly, when applied to 
a system of three qubits, the expression for the moment 
contains all three bipartite concurrence terms, as well as 
the three tangle. In this way this measure captures all 
classes of entanglement, not only bipartite or otherwise. 

Applying the results of Sugita to explain the entan- 
glement behaviour we observe is now straight forward. 
Well below the critical parameter, the ground state Q- 
function will be localized around the solitary fixed point, 
and hence have low degree of entanglement. At the other 
limit, well above the critical parameter, the ground state 
consists of a superposition of two, well-separated (and 
hence almost orthogonal) states. For a bipartite system, 
this implies that the entanglement of formation will be 
1. In between these two limits the distribution is spread 
between the three fixed points, with a greater degree of 
entanglement. As j becomes very large, and the distri- 
butions may be more localized, the point at which the 
ground state distribution is most delocalized will move 
closer towards the bifurcation point. 

In this way, it is the pitchfork structure of the bifurca- 
tion that is vital. The loss of stability of the original fixed 
point, coupled with the emergence of two degenerate sta- 
ble points, as opposed to a single emergent point, which 
is the case for other classes of bifurcation, results in the 
characteristic peak in the ground state entanglement. 

We argue that this result should hold for a general 
bipartite system, and put forth the following conjecture: 

Conjecture 2 (Ground State Entanglement) 

Consider a quantum Hamiltonian H{g) which depends 
smoothly on a parameter g and which acts on a bipartite 
Hilbert space Vi{P) (8) V2{P). The parameter P allows 
one to take a well-defined classical limit. Suppose that 
Ti-ig) is the well-defined classical limit of H{g) and that 
there is a supercritical pitchfork bifurcation of the fixed 
points at the critical parameter, g — gc. Let the von 
Neumann entropy S , be the measure of entanglement. 
Then S{g), the entanglement of the ground state of 
H{g), is a maximum with respect to g at gqc{P) where 
9qc{P) gc in the classical limit. 



IV. SUMMARY 

We have illustrated how the entanglement in the 
ground state of a simple coupled tops model can be as- 
sociated with a bifurcation of the classical fixed points. 
Following this observation, we have argued why this re- 
sult should be generalizable to any quantum system with 
an appropriate classical limit which exhibits a pitchfork 
bifurcation. 

We have already found other instances where our con- 
jecture holds. In particular, for the Dicke model studied 
in Ref. [s^l- This system exhibits a QPT, which corre- 
sponds to the bifurcation. This system is the subject of 
a future article ^35,]. 

The classical bifurcation having a signature - namely 
the entanglement spike - in the quantum regime is not 
that surprising. Most classical characteristics arise in the 
quantum regime, although the signature here - entangle- 
ment - is uniquely quantum. Here we have a correspon- 
dence between the classical and quantum in stationary 
states of the system. Predominantly, quantum-classical 
correspondence has been considered with respect to the 
dynamics, especially in the case of chaotic systems (for 
example Ha, 113, IH, H^, Efll ) • The dynamical generation 
of entanglement is argued to be related to the underlying 
chaos in many-body systems - classical chaos implies a 
greater degree of entanglement. The simple coupled tops 
system offers an excellent test bed for further investiga- 
tion into the relation between chaos and entanglement 
due to its rich (classical) dynamical structure. 
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APPENDIX A: SEMICLASSICAL LIMIT 

To take the appropriate semiclassical limit we first 
consider the correlation function between two operators, 

{X,Y)^{XY)-{X){Y), (Al) 

also known as the covariance. In our case, all operators 
are elements of the SU{2) group of total angular momen- 
tum operators. The scaling of the covariance with respect 
to the total angular momentum eigenvalue, j{j -\- 1), is 

{X,Y)=0{y^JUTT)). (A2) 

i.e., the covariance is of order not exceeding j. Con- 
versely, both (XY) and (X){Y) are 0{j{j -\- 1)). Re- 
expressing Eq. IjAip . and dividing through by j(j + 1) 
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yields 

(XY) 



{X) 



o 



i(i + i) V7(J+i) V7(7+i) VVW+i)/' 

(A3) 

Taking the limit of j oo, meaning + 1) j, 

gives 



{XY) ^ (X) (Y) 



(A4) 



3 3 

. Thus by defining variables as the expectation values 

{Jaot) 



^aot — 



(A5) 



which are simply real numbers, allows the expectation 
values of products of operators to be factorized in the 
scmiclassical limit. The semiclassical dynamics are then 
obtained from the Heisenberg operator equations of mo- 
tion by replacing the operators in the above differential 
equation with these expectation values. Similarly, the 
semiclassical Hamiltonian, E is obtained by taking the 
expectation value of the Hamiltonian, scaling by 1/j, 



E = — — = Lxi + Lx2 + nLziLz2- 
3 



(A6) 
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